Generalized hole-particle transformations and spin reflection positivity in 

multi-orbital systems 
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We propose a scheme combining spin reflection positivity and generalized hole-particle and orbital 
transformations to characterize the symmetry properties of the ground state for some correlated 
electron models on bipartite lattices. In particular, we rigorously determine at half-filling and for 
different regions of the parameter space the spin, orbital and 77 pairing pseudospin of the ground 
state of generalized two-orbital Hubbard models which include the Hund's rule coupling. 



I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tors, a renewed interest has increased in models describ- 
ing strongly correlated electrons. As the simplest one the 
Hubbard model has been intensively studied within the 
analysis of strongly correlated electron systems. Apart 
from it, various generalized Hubbard models have also at- 
tracted considerable interest. "'^ Indeed, by including some 
additional interaction terms into the Hubbard Hamilto- 
nian, either the magnetic ordering phases or the super- 
conducting phases may be stabilized in appropriate re- 
gions of the parameter space. ^ The introduction in these 
models of orbital degrees of freedom'^ enriches the phase 
diagram of the model, because the interplay of orbital 
and spin dynamics can yield a series of novel physical phe- 
nomena, such as unconventional metal-insulator transi- 
tion, novel quantum spin and orbital ordered states, and 
transport property anomalies.* 

Generally, all of these models retain only part of the 
electron-electron interactions which arise as a result of 
the repulsion between two electrons having opposite spins 
and located at the same lattice site. In spite of this crude 
assumption, these models correspond to a many-body 
problem and in general cannot be exactly solved. Ex- 
act solutions are of great importance, since in some cases 
the errors introduced by the approximations may domi- 
nate the results to such an extent that one might end up 
with an incorrect description of the phenomenon under 
study. Considering the exact results with respect to the 
dimensionality of the model, we know that most of them 
have been derived in two limiting case: either in one di- 
mension or in the infinity dimensions. For instance, the 
exact solution of the Hubbard model was given in the one 
dimensional case by means of the Bethe ansatz by Lieb 
and Wu^. The other class of exact solutions belongs to 
the other limiting case, i.e. d=oo, where the dynamical 
mean-field approximation becomes exact. ^ The situation 
gets more complicated as physically interesting lower di- 
mensional cases are considered. d>l rules out the ap- 
plicability of the well-established Bethe ansatz approach, 
while mean-field-like descriptions lead to qualitatively or 
quantitatively incorrect conclusions, because the effects 
of spatial and dynamical fluctuations are not properly 
taken into account. Nevertheless, exact results holding in 



any dimensions are available. We refer to, for instance, 
the Lieb theorem^, the flat-band ferromagnetism® , and 
the Nagaoka theorem.^ However, for extended versions 
of the Hubbard model exact results are still rare.^",^^ 



In this paper we rigorously prove some exact results 
for a generalized Hubbard model. Namely, we consider 
a two-orbital Hubbard Hamiltonian which incorporates 
also the Hund's rule coupling. We notice that the study 
of relevant physical systems requires the use of such kind 
of model whose minimal constituents are the hopping 
term between different, or same, orbitals, and the on- 
site Coulomb and exchange interactions. This is for in- 
stance, to report few of them, the case of non-metallic 
Cu^+compounds^^ (d^ conflguration, one hole in the two 
degenerate eg-orbitals), low spin Ni'^+-^'^(d'' conflgura- 
tion, one electron in eg-orbitals) , as well as Mn'^+ and 
jQj-^g (bigh-spin d* conflguration, one eg electron), 
and in the Ru'*+ ions^^ (d"* configuration, two holes in t2g 
-orbitals). 



Here, our aim is to determine rigorously the symme- 
try properties of the above mentioned Hubbard Hamil- 
tonian by investigating which values of the spin, orbital 
and 77 pairing pseudospin quantum numbers occur in the 
ground state. A physical interpretation and discussion of 
the outcome is also presented. The following conditions 
are assumed for the validity of the results obtained be- 
low: i) the number of electrons is equal to twice number 
of sites in the lattice; ii) the lattice satisfies the connec- 
tivity condition; iii) the hopping amplitude is different 
from zero only for charge transfer between orbitals of the 
same type. 

The scheme of the proof develops into three preliminary 
steps which are based on the use of the property of reflec- 
tion positivity of the Hamiltonian in an opportune range 
of parameters and on the application of continuity ar- 
guments derived by the uniqueness of the ground state. 
Furthermore, the link introduced by some unitary trans- 
formations between the ground states in different regions 
of the space of parameters, allows to recover all the miss- 
ing information about their symmetry properties. 
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II. THE MODEL 

The explicit form of the Hamiltonian for a double or- 
bitally degenerate Hubbard model on a connected bipar- 
tite lattice r is written as: 



where 



H = Ht + Hint, 



Ht = -t 4\ad3\a 
<ij>,X,a 



Hint = Uo'^niX'^nixi + C/ ^ niiani2a+ 



(1) 



(2) 



niii7'ni2o 



Icrdilifdi2gdi2cr 



(3) 



Here dl-)^^ is the creation operator of correlated electrons 
with spin a at site i on orbital A (=1,2), respectively and 
n,A(7 is the number operator for the electrons on the i 
site and A orbital. Moreover, we have used the simplified 
notation a = —a. The parameter t is the hopping ma- 
trix between nearest-neighbor sites, and we assume that 
the A orbitals are not mixed by the hopping; Uo, U and 
J stand for the intra-orbital, intcr-orbital Coulomb and 
Hund's rule exchange interaction, respectively. 

The results presented below are valid for any choice of 
the parameters in the Hamiltonian. Nevertheless, since 
the two orbitals arc equivalent and can be interchanged 
by a properly chosen canonical transformation, we im- 
pose an additional condition on the set of the parame- 
ters, that is Uo = U + J}^'^'^ Hereafter, we will use this 
relation together with the condition that U and J are 
positive. When will be needed, the sign of the Coulomb 
and Hund coupling will appear explicitly. Hence Hint 
becomes, 



Hint — 



(f/ + J) ^ nixiTiixi + n,i„ni2a+ 



where cr arc the Pauli matrices and e{i) = ±1 depending 
to which of the two subparts of the bipartite lattice the 
site i belongs. 

Here S is the usual total spin operator, T is the pseu- 
dospin operator and rj is the so called pairing operator 
introduced by Yang^^, extended to the case of two types 
of electrons. The T operator has properties that are ex- 
actly analogous to the properties of the usual one-half 
spin operator. Indeed, its third component at each site 
assumes the values ^ and — i corresponding to the oc- 
cupied orbitals A = 1 and 2, respectively. Besides, 
takes a fcrmion in the orbital 2, at the lattice site i, and 
transports it to the orbital 1 located at the same lattice 
point. Obviously, Tr corresponds to the opposite pro- 
cess. 

As for S and T, the -q operators generate another SU (2) 
algebra which has as base configurations those with one 
orbital doubly occupied or empty on each site. 

In this model, the total spin operator S commutes with 
the Hamiltonian [S,i?] = and is a good quantum num- 
ber. Another symmetry is that one related to the orbital 
degree of freedom as defined by T. The square of the 
total orbital pseudospin operator T and its third compo- 
nent commute with the Hamiltonian H, since there 
is no hopping between diff'erent orbitals in Ht- On the 
other hand, the square of the total rj pseudospin opera- 
tor and its third component rj^ also commute with the 
Hamiltonian. Hence, these relations imply that these are 
conserved quantities and that the eigenstates of H can 
be labelled as follows: 

!• >= \E, S, Sz,T,Tz,r],riz > ■ 

We also notice that and T+ commutes with H, 
namely the T— multiplets are degenerate in energy. This 
property is not valid for the 77"*" and 77" operators. 

It is important pointing out that these algebras are not 
independent each-other and can be related by means of 
extended hole-particle and orbital type transformations. 

Let us consider the following unitary transformation 
G: 



GdiiiG ^ = dill 



(8) 



Grfjii G ^ = di 



{U - J)^ niiani2a - dl^^diiffdl2^di2a (4) 



Let us introduce the following operators: 

i,(T.<7' A 

i,a,X,X' 



(5) 
(6) 

(7) 



2T 

Gdi2-[G ^ = dill 

Gdi2iG ^ = di2i- 
Under this transformation one has 
GSG"^ =T, 



(9) 
(10) 

(11) 
(12) 



i^a,a' X 



and the Hamiltonian H is transformed as follows: 

GHG-''=G{Ht + Hint)G-^=HG. (13) 
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Here Hq = H(t, J ^ — J), i-e. Hq can be obtained from 
H replacing J with —J, implying that Hq and H are 
related to each other. 

There is also another transformation that will appear 
very useful later for the determination of the symmetry 
properties of the ground state. This is the hole-particle 
transformation of the base operators of creation and an- 
nihilation; 

RdatR-"^ = c^iiT (14) 

RdistR"^ = rfi2T (16) 

Rdi2iR-' = £(i)4i- (17) 
Under this transformation one has 

RSR-^ =T], (18) 
and the Hamiltonian H is transformed as follows: 

RHR-'-=R{Ht + Hint)R-^=HR. (19) 

Here Hr = H(t.J -J,U ~U), i.e. Hr can be 
obtained from H replacing J with — J and U with —U, 
implying that Hr and H are related to each other. 

Our purpose is to determine rigorously the symmetry 
properties of the two-band Hubbard model as defined 
above, by investigating which values of the spin, orbital 
and rj pairing pseudospin quantum numbers occur in the 
ground state. 

As we have pointed out in the Introduction, the proof is 
based on three preliminary steps, followed by the appli- 
cation of the R and G transformations to connect the 
ground states in different regions of the space of parame- 
ters. We will describe how this scheme may allow to get 
rigorous information on the symmetry property of the 
ground state even where the of spin reflection positivity 
is not directly applicable, thus yielding a more general 
description of the ground state symmetry properties. 
Step I: 

Let consider firstly the case derived by applying the 
transformation J ^ — J and U ^ —U on H{U, J) which 
changes H{U,J) H{—U,—J). We remind that this 
transformation is realized by the application of the uni- 
tary transformation R on H{U,J). Having this modi- 
fied Hamiltonian the property of spin reflection positiv- 
ity, it is possible to state^ that the ground state \G)i of 
H{—U, — J), for any positive value of U and J, is unique 
and has total spin S = 0. 

Moreover, for this case it is possible to determine the 
value of the orbital pseudospin quantum number, by 
looking at the relevant configurations which contribute 
to the ground state at values of \U\/t ^ 1. As — f/ is 



large and negative, one ends up with states which con- 
tain only configurations with two electrons occupying the 
same orbital on one site. The total orbital pseudospin for 
this kind of situation is T = 0, since it is identically zero 
on any configuration with two electrons in the same or- 
bital state, i.e. it is a singlet in the orbital space. It is 
worth pointing out that the previous statement can be 
generalized only at the condition that the ground state is 
non-degenerate. This property permits the use of conti- 
nuity like arguments to extend the validity of the result 
from the extreme large negative limit to any finite value 
oi-U. 
Step II: 

Let now analyze the parameter case relative to the Hamil- 
tonian H{U, J) which is the most relevant for systems of 
physical interest. It is worth pointing out that the ground 
state of H{U, J) is unique as it is directly connected to 
that one of H{—U,—J) by the application of the unitary 
transformation R defined in Eqs. 14-19. 
Moreover, by performing a perturbation up to second 
order in U/t, it is possible to show that the low energy 
processes of H{U, J) are described by an eff'ective Heisen- 
berg model for spin one half on an array composed by two 
sublattices A and B defined depending on the sign of the 
magnetic exchange between the spin belonging to the two 
sublattices. The efi'ective Hamiltonian reads as follows: 

Hii = Ji '^{Si(^A)Sj(A) + Si(^B)Sj{B)) + J2 ^ Si(^A)Sj{B) 

where Ji = —J, i.e. the direct Hund coupling and 
J2 = 4i^/([/ + J). For this Hamiltonian, it has been 
proved^'' that the ground state is unique, apart from 
the usual SU{2) degeneracy, and it has total spin S = 
— B\), where A{B) are the total number of lat- 
tice sites belonging to the sublattice A{B), respectively. 
For building the mapping, one has to consider an ar- 
ray where each orbital Aj on the generic j site is associ- 
ated with an effective one, so that — is equal to 

'^\{NxM)-NxAB)) + {NxM)-N^AB))\, Nx,{K) be- 
ing the number of Aj orbitals belonging to the sublattice 
K. 

It is important to stress that it is again the unique- 
ness property of the ground state of Hjj and H{U, J) 
that can permit the use of continuity arguments to state 
that the ground state |G)// of H{U,J) is nondegener- 
ate except for the SU{2) symmetry, and it has total spin 
S^U\A-B\). 
Step III: 

Finally, let analyze the situation obtained by applying 
to the Hamiltonian H{U,J) the transformation R which 
changes J {—J)- By means of this transformation, 
H{U,J) is modified in H{U,—J) where the Hund cou- 
pling assumes the form of antiferromagnetic instead than 
ferromagnetic type interaction. 

As it has been performed in the step II, by means of an 
expansion up to the second order in U/t, the search of 
the symmetry properties of the ground state of H{U, — J) 
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may be carried into the analysis of the ground state for an 
effective Hcisenberg model with two kinds of exchanges. 
The Hamiltonian for the lower energy processes is given 
by the following expression, 

Hiii = Ji '^{Si(^A)Sj(A) + Si(^B)Sj(^B)) + J2 y^'S'j(^)6'j(B) 

»,i id 

where Ji = — J, i.e. the direct Hund coupling with in- 
verted sign and J2 = At^ /{U + J). As for the previous 
case, by using the arguments of Lieb^'^, one can show 
that the ground state is unique except for the SU{2) 
degeneracy and has total spin S = j{\A — B\). Following 
the same procedure of step I, one can state by means 
of continuity arguments that the ground state \G)iii of 
iJ(?7, — J) is unique and has total spin S = \{\A — B\). 
Let remind that the value of A{B) is related to the 
number of orbitals belonging to the sublattice A{B), re- 
spectively. 



(II) an) 




(I) 

FIG. 1. Scheme of the symmetry properties of the ground 
state of the two-band Hubbard model for diflferent regions of 
the parameter space as derived from the procedure of step 
/ - III. 

As a summary in Fig. 1 it has been reported the main 
findings out of the procedure performed above. Each box 

contains the information about the quantum numbers of 
the ground state and to which region of the parameter 
space it refers. The ket \S = a,T =?,r? =?) stands for 
the ground state, whose S,T, and 77 are the values of 
the total spin, orbital and 77 pairing pseudospin quantum 
numbers with a given by ^(lA — 

As one can see in Fig. 1, there are some quantum num- 
bers which cannot be derived from the analysis presented 
in step / — ///, and for this reason the missing values 
have been indicated with a question mark. At this point, 
it is crucial the use of the transformations G and R to 
link the different ground states and the correspondent 
quantum numbers. Indeed, taking advantage of these 
symmetry transformations, one can complete and extend 
the scenario as in Fig. 1 by building up the scheme as it 
is presented in Fig. 2. 



To figure out how to determine the missing informa- 
tion on the ground state, let remind that if we have a 
simultaneous cigenstate of the operator S and of H{U^ J) 
with assigned spin eigenvalue given by s, then the unitary 
operator G transforms this state into an eigenstate of 
H{U, — J) and of T whose pseudospin orbital eigenvalue 
is equal to s too. Then, one can proceed from case II to 
III and extract the information that the ground states for 

([/, J) and H{U, — J) have the same value of the orbital 
pseudospin, i.e. a, so that |G)// = \S = 0,T = a,ri =1 > 
and \G)ni = \S = Q,T = a,r] =? >. 
In the same way, the unitary operator R is mapping 
an eigenstate of the total spin and of the Hamiltonian 
H{U^ J) into an eigenstate with equal total r] pairing 
pseudospin and of H{—U, — J). 

Hence by linking via R the ground states of I and II 
whose one knows the value of the total spin, it is possible 
to deduce the pseudospin values which are now equal to 

and a, respectively. After that, we have full knowl- 
edge of the symmetry character of the ground state in 

1 and II as given by |G)/ = \S = 0,T = 0,r] = a) and 
\G)ii = \S = a,T = a,r] = 0), though it is still missing 
the information on the t] pseudospin in the region HI. 



(II) (III) 




a) (IV) 

FIG. 2. A sketch of the links introduced by R and G be- 
tween the Hamiltonian and their respective ground states as 
discussed in step / — ///. Each box contains now the com- 
plete information about the quantum numbers of the ground 
state in the different regions of parameter space. 

It comes natural at this point to enlarge the range of 
investigation by moving to the case IV which corresponds 
to the region of parameter space where U is negative and 
J is positive as indicated by H{—U,J), where a priori 
we do not have any indication on the symmetry char- 
acter of its ground state. However, the region IV can 
be reached either via I or HI so that we can extract the 
complete information on all the quantum numbers and 
cover the missing points in II and III. The path from I 
to IV turns out to be useful for determining the spin and 
orbital pseudospin values while the way to IV via III can 
yield the information on the rj pairing pseudospin so to 
complete the features for the ground state in these two 
regions, that is \G)iii = \S = 0,T = a,r] = > and 
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\G)iv = |S' = 0, T = 0, ?7 = a >, respectively. 
We would like to point out that there is an alternative 
way to determine the total spin in the ground state which 
does not go through the perturbative arguments used in 
steps II-III and that it only requires the condition of semi- 
positive definiteness of the ground state(s). The theorem 
which refers to this method has been demonstrated and 
extensively discussed in Ref.^^ and it is based on the con- 
dition that a positive definite state and a positive semi- 
definite state built on the same basis vector space have 
a nonzero overlap. The key element of this proof is to 
show that the ground state of the assigned Hamiltonian 
is positive semi-definite and that the ground state of a 
related spin Hamiltonian built on the same lattice and 
on the same Hilbert space is positive definite. Hence the 
condition of zero overlap between those states ensures 
that they have the same value of the total spin. The 
relevance of this way of proceeding is that so doing one 
does not have to prove the non-degeneracy of the ground 
state which is fundamental in the proof that makes use 
of perturbative like arguments. 

For the case of Hamiltonian (1), the proof of the semi- 
positive definiteness of the ground state is not direct be- 
cause one has to use a set of complete and orthonormal 
basis built by means of the operators R and G. Indeed, 
due to the property of spin reflection positivity and to 
the fllling chosen, one can introduce a set of complete 
and orthonormal basis of the form {K|0jj) <8) |?i>^)}, where 

K = R(G) and |(/)^) = Hj.AeQ Axa ^'^^^ ^ being one of 
the possible conflgurations. Then, it is possible to prove 
that states obtained by an expansion on the previous ba- 
sis as \W) = X]a,/3 ^a,/3K|0<!j) ® l?!)^) havc the properties 
of semi-positive definiteness, that is all the eigenvalues of 
the matrix Wa^p are not less than zero. 
At this point it is possible to recover the results of steps 
H-HI by choosing a spin Hamiltonian of the form of i?//- 
Hiii whose the ground state is positive definite and use 
the condition of zero overlap to get the value of the total 
spin. 

Let consider now few comments about the physical con- 
tent of the symmetry character of the ground states in the 
regions I-IV. One important observation is that is 
a novel quantum state which shows coexistence of unsat- 
urated ferromagnetic spin and orbital order on a lattice 
assumed that B,'m the sense that the ground state 
exhibits a value of the total spin and orbital momen- 
tum which is different from zero and scales as the total 
number of lattice sites. Moreover, as discussed in Ref.^^, 
the state \G)i and \G)iv which have r] = a can support 
ODLRO (off diagonal long range order) as it concerns 
the long distances pair correlations. It is striking that 
even in presence of positive Coulomb interaction (i.e re- 
gion IV) one can have an r] pairing in the ground state. 
This can be understood considering that the pairing for 
a multi-orbital system can come in a higher orbital mo- 
mentum state induced by the magnetic exchange so that 
the effect of the local Coulomb repulsion is reduced. 



Another interesting physical property comes from the r] 
pairing symmetry of the ground state of the region II and 
III. Due to the value assumed by rj, it can be shown^°'^^ 
that this state can support coexistence of charge density 
wave (LRO) and superconductivity (ODLRO) in a phase 
which is usually indicated as supersolid. The statement 
is true in the sense that if LRO appears then by symme- 
try the state manifests also ODLRO and viceversa. 
Indeed, it comes that due to the rotational invariancc of 
the ground state in the r] space, the Fourier transform of 
the transverse correlation function at equal time r/+~ (q) 
are proportional to the diagonal one 77^ (q) , thus if one of 
the two has the property to stay finite for some value of 
the momentum q in the thermodynamic limit, the same 
happens for the other one. 



III. CONCLUSIONS 

In conclusion we have shown that the symmetry fea- 
tures of the ground state at half filling of a generalized 
Hubbard model which includes the Hund's rule coupling, 
can be obtained in a large range of the paramater space, 
only by combining the property of spin reflection pos- 
itivity and the use of special unitary transformations. 
In particular it has been possible to extract all the most 
relevant symmetry characteristics related to conservation 
laws of the Hamiltonian that are represented by the or- 
bital, spin and 77 pairing pseudospin operators. We want 
to notice that the use of spin reflection positivity has 
been largely used in literature to rigorously determine 
the symmetry character and the correlation fimctions of 
the ground state for different model Hamiltonians. In 
our work this idea is generalized and enriched in a simple 
scheme which takes advantage of the positivity character 
of the ground state in one region of the parameter space 
and of the use of symmetry transformations to scan a 
larger part of the phase diagram. 

Moreover, the procedure presented above gives also the 
opportimity to extract information on the phase diagram 
of the model in exam, such as the occurrence of a novel 
quantum state with coexistence of spin and orbital or- 
der, the manifestation of ODLRO, and where the ground 
state can be a supersolid. 

Finally, it is worth pointing out that a similar model has 
been exactly solved in the one-dimensional case in the 
case of strong coupling regime, i.e. when the Coulomb 
repulsion is such that the double occupancy of electrons 
on the same site of the same orbital is excluded. 
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